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The Coulomb interaction between the two protons is included in the calculation of 

three-nucleon hadronic and electromagnetic reactions using screening and renormalization 

approach. Calculations are done using integral equations in momentum space. The 

\0 ' reliability of the method is demonstrated. The Coulomb effect on observables is discussed. 
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1. Introduction 



^ ■ The inclusion of the Coulomb interaction in the description of the three-nucleon con- 
^ tinuum is one of the most challenging tasks in theoretical few-body nuclear physics. The 
Coulomb interaction is well known, in contrast to the strong two-nucleon and three- 
nucleon potentials mainly studied in three-nucleon scattering. However, due to its 1/r 
^ ' behavior, the Coulomb interaction does not satisfy the mathematical properties required 
Ch ■ for the formulation of standard scattering theory. There is a long history of theoretical 
work on the solution of the Coulomb problem in three-particle scattering. Some of the re- 
^ ! cent suggestions p!|2|3] have not matured yet into practical applications, while the others, 
^ I based mostly on the configuration-space framework [H[5ll6] . are limited to energies below 
deuteron breakup threshold (DBT). Up to now only few approaches led to the results 
above DBT. Those are configuration-space calculations for proton-deuteron (pd) elastic 
scattering using the Kohn variational principle [7] and the screening and renormalization 
approach in the framework of momentum-space integral equations pi9|10|ll|[T2] : never- 
theless, for the latter method only the present work published in Refs. [TT][T2] uses realistic 
interactions together with fully converged calculations in terms of screening radius and 
two-nucleon and three-nucleon partial waves. 

Section [2] shortly recalls the technical apparatus underlying the calculations and demon- 
strates the reliability of the method. Section [3] gives our summary. 

2. Screening and renormalization approach 

Our treatment of the Coulomb interaction is based on the idea of screening and renor- 
malization proposed in Ref. [T3| for the scattering of two charged particles. The screened 
Coulomb potential of our choice in r-space representation is given by 

«;^(r) = «;(r) e"(^/^)", (1) 
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where w{r) is the proper Coulomb potential, R is the screening radius, and n controls 
the smoothness of the screening. The standard scattering theory is formally applicable 
to the screened Coulomb potential, e.g., the transition matrix is defined via Lippmann- 
Schwinger equation t^(ej + zO) = wji + WRgo{ei + iO)tji{ei + iO) and the corresponding wave 
function is \'4'^^\pi)) = [I + go{ei + iO)t ^{ei + iO)]\pi) , where gQ{ei + iO) is the free resolvent 
and I Pi) is the plane- wave state with momentum pj and energy Cj. Furthermore, as shown 
in Refs. [T3fT^ . the on-shell screened Coulomb transition matrix {pf\tji{ei + iO)\pi) with 
Pf = Pi and wave function diverge in i? — >■ cxd limit, but after renormalization with also 
diverging phase factor ZR^pi) they converge to the proper Coulomb amplitude {pf\tc\Pi) 
and proper Coulomb wave function \ijj'^\pi)), respectively: 

lim {Pfltnici + iO)\Pi)z]^\pi) = (p/|tc|Pi), (2a) 

it— »oo 

hm \^'i^\p.))z-Hp,) = \/^\p.)). (2b) 

The renormalized screened Coulomb amplitude converges to the proper Coulomb ampli- 
tude in general as a distribution. As discussed in Ref . [13] , this is fully sufficient for descrip- 
tion of real experiments and justifies the replacement of lim7j^oo(p/|^i?(ei -l- iO)\pi)z]^^{pi) 
by {pf\tc\pi) in practical calculations [TUfTT] . 

The screening and renormalization approach can be applied to more complicated sys- 
tems, if in the scattering amplitudes the diverging screened Coulomb contributions can 
be isolated in the form of two-body on-shell transition matrix and two-body wave func- 
tion with known renormalization properties ([2D. For the description of (pd) scattering we 
employ Alt-Grassberger-Sandhas (AGS) three-particle scattering equations [15] in mo- 
mentum space 

ulsl^Z) = 6,^G,\Z) + J2'^,.Ti^\Z)Go{Z)Ui^Hz), (3a) 
uiS'iZ) = G,\Z) +Y.n^\Z)G,{Z)Ui^^{Z), (3b) 

cr 

with 6pa = 1 — ^i3a, Gq{Z) being the free resolvent, Ta^\z) the two-particle transition 
matrix derived from nuclear plus screened Coulomb potentials, and uj^^{Z) and U^^\z) 
the three-particle transition operators for elastic/rearrangement and breakup scattering; 
their dependence on the screening radius R is notationally indicated. As demonstrated in 
Refs. [5|11II12] . the three-particle transition operators can be decomposed into long-range 
and Coulomb-distorted short-range parts 

U^0aiZ) = Sp^T^^iZ) + [1 + T^^{Z)Gf\z)]uj,^\z)[l + Gi''\Z)T^^{Z)], (4a) 
U^a\z) = [1 + V(Z)Go(^)]C^(Z)[l + Gi'^\Z)T^^{Z)] (4b) 

with the channel resolvent G^J^\z), proton-proton (pp) screened Coulomb transition ma- 
trix TpR{Z), two-body transition matrix T^^{Z) derived from the screened Coulomb po- 
tential between spectator and the center of mass (cm.) of the remaining pair, and with 
the reduced short-range operators U^J{Z) and UqJ{Z). On-shell matrix elements of the 
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operators (jlj) between two- and- three-body channel states |0a(qj)t'ai) and |0o(p/q/)t'Oj) 
with discrete quantum numbers u^., Jacobi momenta Pj and q^, and energy Ei, do not 
have a i? — >■ cx) hmit. However, the quantities diverging in that hmit are already isolated 
in Eqs. @ and are of two-body nature, i.e., the on-shell T^^{Z) and the initial/final state 
screened Coulomb wave functions. Those quantities, renormalized according to Eq. ([2]), in 
the R ^ oo limit converge to the two-body Coulomb scattering amplitude and to the 
corresponding Coulomb wave functions, respectively, thereby yielding the pd scattering 
amplitudes in the proper Coulomb limit 

{(f)(3{cif)iyi3f\Upa\(f)a{qi)jya,) = 5/3^ (0a (q/) Z^a^ I (^i ) '^"^ ) 

+ il^J^PR iM^f^Pf I [U^pa (E^ + ^0) (5a) 
- 6p^T^^iE, + ^O)]|0,(q,>„jzJ(gO}, 
(0o(P/q/)z^O;|f/oa|0a(qi)z^a,) = lim {2;^^(p/)(0o(p/q/)z/oy.| 

xf/(f)(E, + zO)|0„(q,)z/.,)2j(g.)}, 

where the relation between Eqs. ([5]) and (jlj) is given in Refs. [TTIT2] . The renormalization 
factors ZaR and zr are diverging phase factors defined in Ref. [T3llllfl2j . The R ^ oo 
limit in Eqs. ([5]) has to be calculated numerically, but due to the short-range nature of the 
corresponding operators it can be reached with sufficient accuracy at rather modest R if 
the form of the screened Coulomb potential ([T]), in particular the parameter n controlling 
the smoothness of the screening, has been chosen successfully. We prefer to work with 
a sharper screening than the Yukawa screening {n = 1) of Refs. piTO] . We want to 
ensure that the screened Coulomb potential WR^r) approximates well the true Coulomb 
one w{r) for distances r < R and simultaneously vanishes rapidly for r > R, providing a 
comparatively fast convergence of the partial- wave expansion. In contrast, the sharp cutoff 
(n — *■ cxo) yields an unpleasant oscillatory behavior in the momentum-space representation, 
leading to convergence problems. In Refs. [TTIT2] we found the values 3 < n < 6 to provide 
a sufficiently smooth, but at the same time a sufficiently rapid screening around r = R; 
n = 4 is our standard choice. The screening radius R sufficient for convergence in Eqs. ([5]) 
is considerably larger than the range of the strong interaction. As a consequence, the 
calculation of the three-particle transition operators for nuclear plus screened Coulomb 
potentials requires the inclusion of partial waves with angular momentum much higher 
than required for the hadronic potential alone. This problem can be solved in efficient 
and reliable way either by using the perturbative approach for high two-particle partial 
waves, developed in Ref. [TB], or even without it as discussed in Ref. [T7|. More details on 
the practical implementation of the screening and renormalization approach as well as the 
extension to three-nucleon electromagnetic (e.m.) reactions are presented in Refs. [TT][T2] . 

The internal criterion for the reliability of our method is the convergence of the observ- 
ables with screening radius R employed to calculate the Coulomb-distorted short-range 
part of the amplitudes in Eqs. Figures [T] and [2] show characteristic examples for pd 
elastic scattering and breakup. The hadronic interaction is the realistic coupled-channel 
potential CD Bonn -f- A, allowing for single virtual A-isobar excitation [18]. In most 
cases the convergence is impressively fast; the screening radius i? = 20 fm is sufficient. 
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Figure 1. Convergence of the pd elastic scattering observables with screening radius R. 
The differential cross section, proton analyzing power Ay{N) and deuteron analyzing 
power T21 at 3 MeV proton lab energy are shown as functions of the cm. scattering 
angle. Results obtained with screening radius i? = 10 fm (dotted curves), 20 fm (dash- 
dotted curves), and 30 fm (solid curves) are compared. Results without Coulomb (dashed 
curves) are given as reference for the size of the Coulomb effect. 
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Figure 2. Convergence of the pd breakup observables with screening radius R. The 
differential cross section in selected kinematical configurations at 13 MeV proton lab 
energy is shown as function of the arclength S along the kinematical curve. Notation of 
curves as in Fig. [H 



The exceptions requiring larger screening radii are the pd elastic scattering observables 
at very low energies and the breakup differential cross section in kinematical situations 
characterized by very low pp relative energy Epp, i.e., close to the pp final-state interaction 
(pp-FSI) regime, as shown in Fig. [31 In there, the pp repulsion is responsible for decreasing 
the cross section, converting the pp-FSl peak obtained in the absence of Coulomb into a 
minimum with zero cross section at pf = 0, i.e., for Epp = 0. Such a behavior is seen 
in the experimental data as well p^][20] . The slow convergence close to pp-FSl is not 
surprising, since the renormalization factor Zfi{pf = 0) itself is ill-defined, indicating that 
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Figure 3. Convergence of the pd breakup observables with screening radius R. The 
differential cross section for pd breakup at 13 MeV proton lab energy in the pp-FSl 
configuration is shown as function the relative pp energy Epp. Results obtained with 
screening radius i? = 10 fm (dotted curve), 20 fm (dashed-double-dotted curve), 30 fm 
(dashed-dotted curve), 40 fm (double-dashed-dotted curve), 60 fm (solid curve), and 
results without Coulomb (dashed curve) are compared. 

the screening and renormalization procedure cannot be applied at pj = 0. Therefore an 
extrapolation has to be used to calculate the observables at pf = 0, which works pretty 
well since the observables vary smoothly with pf. 

Furthermore, Ref. |21] makes a detailed comparison between the results for pd elastic 
scattering obtained by the present technique and those of Ref. obtained from the 
variational solution of the three-nucleon Schrodinger equation in configuration space with 
the inclusion of an unscreened Coulomb potential between the protons and imposing 
the proper Coulomb boundary conditions explicitly. The agreement, across the board, 
between the results derived from two entirely different methods, clearly indicates that 
both techniques for including the Coulomb interaction are reliable. 

3. Summary 

We have shown how the Coulomb interaction between the protons can be included into 
the momentum-space description of proton-deuteron scattering using old idea of screening 
and renormalization [13]. The theoretical framework is the AGS integral equation [T5] . 
The calculations are done on the same level of accuracy and sophistication as for the 
neutron-deuteron scattering. Our practical realization differs enormously from the one 
of Refs. pun], even the form of three-particle equations including screened Coulomb is 
different. We use modern hadronic interactions whereas the calculations of Refs. PITO] 
were based on quasiparticle equations with rank-1 separable potentials and in addition 
approximated the screened Coulomb transition matrix by the screened Coulomb potential. 

Compared to the configuration-space treatment [7], the results for the elastic pd scat- 
tering agree very well over a wide energy range [21]. Although at very low energies the 
coordinate-space methods remain favored, at higher energies and especially for three-body 
breakup reactions our momentum- space treatment is more efficient, and so far the only 
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one to show first results for realistic interactions |22j . 

A realm of applications to the rich pd data base for elastic scattering and breakup may 
be seen in Refs. [TT|12lll9ll^ . The conclusion is that in elastic pd scattering the Coulomb 
effect is important at low energies for all kinematic regimes, but gets confined to the 
forward direction at higher energies. In pd breakup and in three-body e.m. disintegration 
of ^He the Coulomb effect is extremely important in kinematical regimes close to pp-FSl. 
There the pp repulsion converts the pp-FSl peak obtained in the absence of Coulomb 
into a minimum with zero cross section [20] • This significant change of the cross section 
behavior has important consequences in nearby configurations where one may observe 
instead an increase of the cross section due to Coulomb [19]. However, some of the long- 
standing discrepancies between experiment and theory like the space star anomaly in pd 
breakup are not resolved by the inclusion of the Coulomb interaction. 

Finally, the screening and renormalization approach for including the Coulomb inter- 
action is extended to four-nucleon scattering [23] . 
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